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5— ( ' Wc study the photoproduction of an isolated photon and a jet based on a code 

of partonic event generator type which includes the full set of next-to-leading order 
corrections. We compare our rcsuhs to a recent ZEUS analysis in which an effective fc^ 
of the incoming partons has been determined. Wc find that no additional intrinsic kT 
is needed to describe the data. 
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1 Introduction 

Inclusive photoproduction of an isolated prompt photon has been measured at HERA |j|] and 
compared to theoretical predictions |2|, ^, |^ recently. The photoproduction of an isolated 
prompt photon and a jet has been studied less extensively [^, and only partial next-to- 
leading order codes were available up to now [^, ^. However, the presence of the jet allows 
the definition of observables which may provide a more detailed picture of the underlying 
partonic processes than in the inclusive case. In particular, it allows to extract information 
on the transverse momentum kx of the initial state partons. 

The introduction of an effective initial state transverse momentum (fcy) [^ has become a 



common practice in recent years p, 0, 10 1 and has served to bridge large gaps between certain 
data and NLO theory in fixed target and hadronic collision experiments of prompt photon 
production. It has been argued that, in addition to the "intrinsic" transverse momentum 
of the initial state partons of a few hundred MeV due to the finite size of the proton, 
initial state soft gluon radiation can generate sizeable transverse components of the parton 
momenta. Considerable progress has been made in recent times to account for multiple 



soft gluon emission effects in resummation formalisms [11[-|14], but the phenomenological 
application of such calculations has not reached a mature state yet. Therefore, simplified 
phenomenological models for kx effects have been employed which introduce a Gaussian 
smearing of the parton transverse momenta, the mean value of this effective kx being fitted 
from data. However, these procedures are rather ad hoc, and various experiments lead to 
quite different conclusions ||l^, |l6[ on the necessity to introduce an effective kx- 

Recently, the ZEUS collaboration reported on the analysis of prompt photon plus jet 
photoproduction data |17|. From the comparison of the data to results obtained with the 



Monte Carlo program PYTHIA |18], it has been concluded in |17| that a total effective {kx) 
value of about 1.7 GeV is necessary to fit the data. This value is composed of the mean value 
of the intrinsic parton momentum in the proton, (A;™*''), and the parton shower contribution 
to {kx)- The two contributions are combined assuming that the overall distribution is 
Gaussian. The parton shower contribution was found to be approximately 1.4 GeV, and for 
(fc™*'') a value of (A;™*'') ~ 1.25 GeV has been fitted. However, it was found that the same 
data are already well described by HERWIG [ p^ ] with a default parton intrinsic transverse 
momentum of zero. 

In this paper, we compare NLO QCD results to the data presented in |l^] and discuss 
in detail possible sources of differences between theory and data. In particular, we address 
the question whether an extra effective transverse momentum {kx) - in addition to the one 
already contained in fixed order NLO QCD - is necessary to fit the ZEUS data. We found 
that our NLO QCD prediction for the quantities which served to determine {kx) in [|l^] 
describes the data without taking into account any extra kx effects. We recall that the 
same is true for a large set of prompt photon data in hadronic collisions and fixed target 
experiments [^, 16, 20|. 



Our results have been obtained with a program of partonic event generator type which 
includes the full next-to-leading order corrections to all the partonic subprocesses, as well 
as the box contribution 7(7 -^ ^g. The program already has been used to study inclusive 
prompt photon photoproduction Q , and has been extended to allow also the calculation of 



prompt photon plus jet photoproduction in the present work. 

The paper is organized as follows. In section 2, we first describe the main lines of 
the calculation our code is based on. Then we discuss the effect of certain kinematic cuts 
and infrared sensitive limits which will be relevant for the comparison to the ZEUS data. 
In section 3 we give numerical results for the isolated photon plus jet cross section and 
compare to ZEUS data. Section 4 contains the study of /cr-sensitive observables. We 
show our results for the p± and A</> distributions used by ZEUS to determine {kx) and 
discuss their implications. An appendix contains details of the discussion of symmetric cuts 
contained in section 2. 

2 Description of the method 

2.1 General setting 

As the general framework of the calculation already has been described in detail in [^, we 
will give only a brief overview on the method here. 

In photoproduction events, the electron acts like a source of quasi-real photons whose 
spectrum can be described by the Weizsacker- Williams formula 

r(y) = ^|l±il^ll!ln^I^4^-^ii^| . (1) 

^ 27r [ y miy^ y J 

The quasi-real photon then either takes part directly in the hard scattering process, or it 
acts as a composite object, being a source of partons which take part in the hard subprocess. 
The latter mechanism is referred to as resolved process and is parametrized by the photon 
structure functions Fai^{x^,Q'^). Thus the distribution of partons in the electron is a 
convolution 

Fa/e{xe, M) = dy dx^ f^{y) F^/^{x^, M) 6{x^y - Xe) (2) 

where in the "direct" case F^^/^^x^, M) = da-yS^l — x^). 

Similarly, a high-p^ photon in the final state can either originate directly from the hard 
scattering process or it can be produced by the fragmentation of a hard parton emerging 
from the hard scattering process. Thus we distinguish four categories of production mech- 
anisms, as illustrated in fig.|l|: 1. direct direct 2. direct fragmentation 3. resolved direct 
4. resolved fragmentation. We implemented the full set of next-to-leading order corrections 
to all four subprocesses. Note that the photon structure functions behave for large Q^ as 
Fa/.y{x,Q'^) ~ InQ^/A^ ~ l/as{Q'^). An analogous argument holds for the fragmentation 
functions, such that a consistent NLO calculation requires the inclusion of the full 0{as) 
corrections not only to the "direct direct" part but also to the resolved and/or fragmen- 
tation parts^ The matrix elements for these NLO corrections have been taken from the 
literature pl|, |22|, 23]. We also included the box contribution gj -^ (77, which is formally 



NNLO, but known to be sizeable [|4|, |2§. 



^In the following, we will denote by "resolved" the sum of resolved direct and resolved fragmentation 
parts and by "direct" the sum of direct direct and direct fragmentation parts. 
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Figure 1: Examples of contributing subprocesses at leading order 



It is well known that at NLO the attribution of certain finite terms to either of the four 
categories listed above is scheme dependent and therefore only the sum of all parts has a 
physical meaning. 

In order to isolate the infrared singularities appearing in the calculation at next-to- 
leading order, a modified phase-space slicing method has been used Q, [20[| . 

For the definition of jets, the fc^-algorithm p^] has been employed, as has been done 
by the ZEUS collaboration. 

To single out the prompt photon events from the background of secondary photons 
produced by the decays of light mesons, isolation cuts have to be imposed on the photon 
signals in the experiment. We use the following isolation criterion: In a cone C centered 
around the photon direction in the rapidity and azimuthal angle plane, defined by 



C = {(r/, </>)! (r? - r/^)' + (</>- 07)' <R^] . 



the amount of deposited hadronic transverse energy EjP'^^R) is required to be smaller than 
a fraction €c of the photon transverse momentum. We use ec = 0.1 and R = 1 to match the 
cuts of the ZEUS collaboration. 

Since isolation cuts impose additional phase space restrictions, the question has been 
raised if factorization properties valid for the inclusive case might be spoiled by isolation. 
This issue is discussed in detail in [31, 32 1. Here we just note that the above criterion, being 



based on boost-invariant transverse energies, leads to infrared safe cross sections. 

2.2 Effects of kinematic cuts, infrared sensitive limits 

As a consequence of phase space restrictions induced for example by kinematic cuts, certain 
observables calculated at fixed order may show instabilities (integrable singularities) at some 
critical point of phase space |3C, 31|. Examples will be discussed below. 

Symmetric cuts on Ej^ , Ej, 

In order to restrict both prompt photons and jets to well-measured kinematic regions, the 
ZEUS collaboration has required Ej. > 5GeV and Eij^ > 5GeV in their analysis [ p^ . 



However, as it is well known ||3^, 34 1, symmetric cuts on the transverse energies should 
be avoided since they select a region where fixed order NLO QCD loses its predictive 
power. This can be seen by considering the total cross section as a function of Ej^^^^ 
as shown in fig.|2|. The cut on the transverse energy of the photon has been fixed to 



Ej, j^jj^ = 5 GeV. The figure shows that the NLO QCD prediction varies rapidly around the 

^T.min ~ -^T,min- 



point Ei^ ■ = EZ ■ In particular, it is striking that the cross section is not immediately 



monotonically decreasing for increasing E^^^^. The reason for the dip at Ej^^^-^^ = Ej.^^^^ 
is that at this point the phase space for the unobserved parton in the real corrections is 
severely restricted, such that the (negative) virtual corrections are not sufficiently balanced 
by the (positive) real correctionsQ As a result we observe a violation of unitarity in the range 
5.0 GeV< E!^^-^^ < 5.25 GeV, where the cross section increases when the available phase 

^We checked by varying the phase space shcing parameter that this effect is not an artifact of the slicing 
method. 



space decreases. This phenomenon, that we discuss in detail in the Appendix, is an artifact 
of the fixed order calculation. In this particular region around Ej^^:^^ = Ej,^-^^, the higher 
order corrections should be resummed to produce a curve that should look like the dashed 
curve in fig.||. This resummation is not considered in this paper (but the large logarithms 
responsible for this behaviour are identified, see Appendix A) because the unphysical effect 
at -Ejfj^jj^ = Ej, ^jj^ is in general numerically weak. However, we shall explore the effect of 
different choices of E^^^^ when calculating the cross section corresponding to ZEUS data. 
It is clear that an asymmetric cut in the experiment would avoid any trouble. 
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Figure 2: Nonisolated total 7+ jet cross section for Ej, > 5 GeV as a function of E^^^^. The 
photon energy range is 0.2 < y < 0.9, the rapidities are integrated over —1.5 < rj^^* < 1.8 
and -0.7 < r]t < 0.9. 



Definition of x^, limit x^ — > 1 

Care also has to be taken in the treatment of the variable x^ near the boundary x^ ^ 1. 

At leading order, the fraction x^ of the incoming photon longitudinal momentum is defined 

as 

_ ET,{e-^^+e-^^) 

^^ ~ 2E^ ^^^ 

where r/3 and 7^4 are the rapidities of the two outgoing partons and E^ = yEe- At next-to- 
leading order, the "true" x-y for 2 ^ 3 processes is given by x^ = J2i=3{^Tie~^^)/{2Ey), 
but since particle 5 is not observed, x^ is no longer an experimentally accessible observable. 
Hence, to obtain an estimate of the longitudinal momentum fraction of the incoming photon, 
the variable x"!'^, defined as 



7 



xf = ""^^ tj"^ ' (4) 



lEy 



is frequently used. However, the above definition of x!? may lead to infrared sensitive cross 
sections ||3^. Fixing E'^ and E^ strongly constrains the phase space of parton 5 when 
^obs gQgg ^Q Qj^Q This is reflected by log(l — x°^^) terms generated by the higher order 
corrections. In general, x'^'^ is integrated over, typically in the range 0.75 < x°^'' < 1, 
such that the singular behaviour is sufficiently smoothed out for a fixed order calculation 
to remain valid. However, if the width of the last bin at x°^^ — > 1 is small, the shape of 
the theoretical prediction near x° = 1 is considerably dependent on the bin width. For 
this reason, the variable x^^ = Ej,{e~^'' + e"'''^" )/{2Ey), which has a smoother behavior 
for Xy -^ 1, has been proposed in ||3^. It also has the advantage that it does not depend 
on Ei!j^ and thus does not contain an uncertainty from the jet energy reconstruction. Note 
that x^^ can be bigger than one, as illustrated in fig.^. 

For the analysis presented in |17], the ZEUS collaboration did not use x°^^, but the 
variable x'^^"'^, defined as 



^ 2E,yJB 2E,yJB 






where the sum is over the prompt photon and the contents of the jet, and y is the energy 
fraction of the incoming photon, reconstructed with the Jacquet-Blondel method: 

Here the sum is over all energy- flow objects in the event, each of which is treated as due 
to a massless particle with energy E and longitudinal momentum component Pz- From 
an experimental point of view, x™^"** is convenient since a lot of calorimeter calibration 
systematics cancel out. The difference between definitions (0) and (0) for x^ depends on 
the definitions of Ei^ and rp^^ in (Q) . If the Snowmass conventions 

^T -1^ ^Ti , r - — =7i^ , r - — =7^^ (7) 
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Figure 3: Behaviour of da/dx^^, da / dx'T''^"-'^ and da/dx°^'^ for the isolated 7 + jet cross 
section. 



are used, one has 
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For our analysis, where maximally two niassless partons form a jet, the difference between 
x'^'^ and x^^""^ is numerically indistinguishable (see fig. S). 

Obviously, what has been said above on the limit x^ ^ 1 is valid for x. 
well. 



meas ^ -i „„ 



3 Numerical results and comparison to ZEUS data 

For the numerical studies we use the rapidity cuts —1.5 < r?-'^* < 1.8 and —0.7 < rf < 0.9 
chosen by ZEUS to restrict the photon and the jet to weh measured kinematic regions. We 
also use the cuts Ei^^:^^^ = Ej.^^^ = 5 GeV in order to match the ZEUS cuts, although such 
symmetric cuts select a region where fixed order perturbative QCD shows instabilities, as 
discussed in detail in section 2^ and in the appendix. From fig. |2| one can estimate that 



the uncertainty in the theoretical prediction introduced by using these cuts is of the order 
of 6-10% for the total cross section. 

We set Qmax ~ 1 GeV^ for the virtuality of the photons in (||) and restrict the photon 
energy to the range 0.2 < y < 0.7 as reconstructed by ZEUS. We take the MRST2 |p7|] 



parametrization for the parton distributions in the proton, AFG |25] for the parton dis- 



tributions in the photon, BEG |29] for the fragmentation functions into photons, nj = 4 
flavors, and for as(/i) we use an exact solution of the two-loop renormalization group equa- 
tion, and not an expansion in log(/x/A). Unless stated otherwise, the scale choice fi = pj, 
has been adopted and M and Mp have been set equal to /i. The rapidities refer to the ep 
laboratory frame, with the HERA convention that the proton is moving towards positive 
rapidity. 

Fig. ^ shows the relative magnitude of the NLO result to the leading ordei]^ result for the 
isolated photon rapidity distribution. As the box contribution has the same kinematics as 
the leading order direct direct part, we included the box contribution in this part, such that 
"direct direct" always means direct direct plus box. Since the higher order corrections to the 
direct direct term are negative and those to the resolved part are small (for the scale choice 
/i = pj,, see fig.|5|), the overall NLO result is lower than the leading order result. We note 
that the higher order corrections to the direct part are most negative for E^^^^ = Ej,^^^, 
for asymmetric cuts the absolute value of the negative contributions decreases. 

In order to single out the events stemming from a direct photon in the initial state, the 
cut x'T'^"'^ > 0.9 has been applied by ZEUS. However, it has to be kept in mind that beyond 
leading order, the direct part has a non- negligible tail away from x™^"* = 1, and that the 
contributions from the resolved part to the bin x^^""^ > 0.9 are rather large, as can be seen 
from Fig. ^(a). Fig.^(b) compares the sum direct + resolved with the cut x'T'^"'^ > 0.9 (0.85) 
to direct direct only with no cut on x^^""^ for the photon rapidity distribution. One observes 
that there is a significant difference in the shape between the direct part and the result 
obtained by imposing the cut x™^"** > 0.9 to the full set of subprocesses, the latter being 
higher in the backward region and lower in the forward region. Hence imposing a stringent 
cut on x^^""^ cannot be considered as fully equivalent to singling out the direct events only. 



^By "leading order" we mean the sum of the lowest order terms of all four subprocesses direct direct, 
direct fragmentation, resolved direct and resolved fragmentation. 
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Figure 4: Magnitude of leading order and NLO results for the isolated cross section 
da^'^^^^/drj"' . Isolation with €c = 0.1, R = 1, jet rapidities integrated over —1.5 < rj^^* < 1.8 . 
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Figure 5: Magnitude of LO and NLO contributions for resolved and direct parts separately, 
with isolation (cc = 0.1, R = 1). 
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3.1 Comparison of photon rapidity and x^ distributions to ZEUS data 

Fig. shows the distributions da/dx^^"''^ and the photon rapidity distribution for various 
scale choices together with the ZEUS data. One can see that the theoretical uncertainty due 
to overall scale variations, i.e. varying /i between Py/2 and 2pj, and setting Mp = M = /i, 
is very small. For the distribution in x™*^"* one observes that theory overshoots the data 
substantially in the last bin, whereas the data are described well in the remaining bins. In 
da/drf there is a rather large discrepancy between theory and data in the forward region. 
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and // = 2pj, (always Mp = M = fi) together with ZEUS data. The errors on the data are 
statistical only. 

Reasons for these discrepancies may be the following: 
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1. The data are at the detector level, i.e. they are not fully corrected for energy losses in 
the detector and other effects. Therefore data corrected to the hadron level may differ 
in a non-negligible way from the ones displayed in fig^. Moreover our calculations 
are performed at the parton level and hadronization effects may also lead to sizeable 
corrections at these low values oi p!^ (see also point 4. below). 

2. There is an uncertainty in the reconstruction of the "true" photon energy y with the 
Jacquet-Blondel method (cf. eq. (|6|)). In general, corrections for detector effects and 
energy calibration tend to shift the range of y to slightly higher values as compared 
to the original range of y^^ ■ To study the effect of an uncertainty in y, we varied 
the cuts on y, increasing the lower cut from ymin = 0.2 to ymin = 0.25, respectively 
increasing the upper cut from ymax = 0.7 to ymax = 0.8, as shown in fig. ^. 

3. Concerning da / dx'T'^'^^ : 

There is an uncertainty in the energy and pz values entering in the definition (^) 
of x^*^"*. Taking into account this uncertainty in the calculation would introduce a 
Gaussian smearing in E and pz which would smooth out the theoretical curve near 

meas 1 

x^ ^ -L . 

4. As already discussed in [Q, a non-negligible part of the hadronic energy measured in 
the isolation cone may be due to underlying events, especially in the forward region. 
Therefore it is possible that experimentally, the isolation cut removes events which 
pass in the theoretical (parton level) simulation since there the underlying event con- 
tamination is not present. It has been estimated [Q] that this can have an effect of 
the order of 30% in the forward region. 

Fig.^ again demonstrates (see also fig.§) that the cross section is very sensitive to a 
variation of E^r^^:^^ around the value Ej^^:^^ = E^.^^^. It can be clearly seen that increasing 

the cut on Ei^ from 5 GeV to 5.5 GeV does not lead to a decrease in the cross section, as 



explained in section 2.2. For the rapidity distribution the theoretical uncertainty introduced 



by using symmetric cuts can be as large as 30%. 

An intrinsic (kx) would not improve the agreement between theory and data for the 
photon rapidity distribution, since it would shift the theoretical curve upwards without 
changing its shape. We note also that the data for da/drf and da jdx^^"''' are in general 
well described by PYTHIA 6.129 with default (/c^^) values |]l^]. Only in the backward 
region in the photon rapidity distribution, PYTHIA is lower than the data, the same trend 
is visible in the NLO theory curve. 



Figure IC shows the contributions from direct and resolved initial state photons sepa- 
rately. As the those parts separately are highly scale dependent, results for two ("extreme") 
scale choices are shown. The figure for the rapidity distribution illustrates that the state- 
ment which part - i.e. resolved or direct part - dominates in the forward region, is scale 
dependent. 
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Figure 8: da"'^^^^ /dx'T'^'^^ and photon rapidity distribution da^'^^'^^ /dif for different cuts 
on the photon energy: 0.2 < y < 0.7 (default), 0.25 <y < 0.7, 0.2 < y < 0.8. 
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Figure 10: Contributions from resolved and direct initial state photons separately for two 
different scale choices. 
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4 Results for kr - sensitive observables and comparison to 
ZEUS data 

In this section we compare our results with the ZEUS data for the kx - sensitive observables 
p±_ and A0 which served to determine {kx) in |17]. The quantity p±_ is the momentum 



component of the photon perpendicular to the jet direction, defined as 

Pl = \p^t^P^''\/p^t' (8) 

Here p^ is not the Snowmass Ej^ , but the modulus of the transverse components of the 
vector sum of the particle momenta which form the jet. Ac/) is the azimuthal acollinearity 
between the photon and the jet. To define A^, the Snowmass definition (see eq. (|^) for 
the jet azimuthal angle has been used by ZEUS. In the experiment, a cut of x"!^^""^ > 0.9 
has been imposed in order to select events which are predominantly from direct incoming 
photons, thereby suppressing eventual contributions to (kx) from the resolved photon. The 
data for these quantities have been corrected to hadron level. As the value of (fcr) affects 
mainly the shape of the distributions, the cross sections have been normalized to one in 
order to minimize uncertainties in the calorimeter energy scale. 



Fig. |ll| shows the p± and Acf) distributions together with the ZEUS data. One can see 
that NLO QCD does describe the ZEUS data well without any intrinsic kx or resummation 
taken into account. A reason for the fact that no need for resummation near A^ = vr or 
p_i_ = is visible might be that the bin size is relatively large. 

Since there is an uncertainty in the reconstruction of x!^^""^ , we investigated the effect 
of shifting the cut on x'T''^"'^, as illustrated in fig.|l^. The reason why the result for the last 
bin p_L — > respectively A(j) -^ vr decreases if the lower cut on x™*^°* is decreased is related 
to the fact that for x™^'*^ > 0.8 more inelastic processes contribute than for x"!^^""^ > 0.9, 
such that the ratio of nearly back-to-back events to the full cross section decreases. 

Note that we used symmetric cuts E^^^^ = Ej, ^:^^^ = 5 GeV in order to match the 
ZEUS cuts. As discussed above, the theoretical result shows instabilities for symmetric cut 
values. Therefore we investigated the effect of varying Ej^^:^^ by 10% around the critical 
point -Ej'rnin = 5 GeV, as shown in fig. 12. One can see that the result varies by less than 



10% in the normalized cross sections. The effect of an uncertainty in the cuts on the photon 
energy y has been found to be negligible in the area normalized plots. 

Hence we can say that the agreement between fixed order NLO QCD and data is very 
good in the case studied here. No kx in addition to the one already implicitly accounted 
for by including the full set of NLO corrections is needed. We recall that he same was true 
for the inclusive case Q, but of course there the sensitivity to kx is smaller. 

Further we note that the intrinsic kx (i.e. the component of {kx) that does not include 
the parton shower contribution) determined by the ZEUS collaboration to fit best the data 
has been obtained by comparing to simulations with PYTHIA 6.129. On the other hand, 
it was also found that the data are already well described by HERWIG 6.1 with the default 
parton (fc^*'') of zero |l^. The reason might stem from differences in the way the parton 
showering is implemented in HERWIG and PYTHIA. HERWIG does not use a sharp lower 
cut-off in the shower evolution below which PYTHIA relies on a suitable value of (/c™*''). 
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Figure 11: Normalized cross section differential in p± and A0 for x!"^"* > 0.8,0.85,0.9. 
The inner error bars represent the statistical errors only, the outer bars represent the total 
uncertainty. 
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Figure 12: Normalized cross section differential in p± and Acp for Ej^^^^^ = 4.5, 5 and 
5.5 GeV with x!""^"^ > 0.9 . The inner and outer error bars correspond to statistical and 
total errors, respectively. 
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5 Conclusions 

We have presented a full next-to-leading order calculation for the photoproduction of an 
isolated photon accompanied by a jet, based on a program of partonic event generator type. 

We found that the variation of our result due to scale changes is very weak. We also 
investigated the sensitivity of the cross section to variations of the photon energy, to different 
cuts on Xry and to the cut E^r^^:^^ on the jet transverse energy. In particular, we discussed 
in detail the instabilities arising at fixed order NLO QCD if symmetric cuts on the photon 
and jet transverse energies are chosen. 

In comparing the x^ distribution to ZEUS data we found that the agreement is good 
except near x^ — > 1. For the photon rapidity distribution we found that our prediction is 
below the data in the backward region and overpredicts the data in the forward region. We 
discussed in detail possible sources of discrepancies between theory and data. 

We also compared our results to recent ZEUS data for fc-r - sensitive observables. Since 
the comparison has been done for normalized cross sections, the effect of variations in the 
photon energy, x^"* or Ej^^:^^ turned out to be very small. We found that NLO QCD 
already describes these data well without introduction of any additional intrinsic kx- 

Our results are in contrast to the results found in the E706 fixed target experiment p|, 
where the conclusion was that high-p^^ prompt photon and vr'' data could not be well de- 
scribed by NLO perturbative QCD without the incorporation of substantial extra kx effects. 
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A Appendix 

In this appendix we identify the origin of the pathological behaviour of the cross section 
in the vicinity of Ej,^:^^ for symmetric cuts on the photon and jet transverse energies 
(see fig. 0). Our aim is to isolate only the leading behaviour, i.e. we work in the double 
logarithmic approximation and drop subleading contributions. 

We consider the subprocess pi+P2 — > P7 +P4 +P5 where ps is the momentum of a gluon. 
This gluon represents a NLO correction to the Born subprocess pi + P2 ^ P7 + Pa and in 
the soft approximation {\p^\ <C \P'y\) this NLO correction can be written as 

CT^^'^ ~ j dcP^dpT, PtI-'' In (^^) ^B°™(pi +P2^Py+ PA) . (9) 

The cut Ej^^:^^ < E!j^ (= pta in the approximation we consider here) restricts the gluon 
phase space to values of pT5 smaller than pT-y and depending on (j)^: p™|^ = PtsI'/'s)- ^^^ 
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The 1/e^ term is cancelled by the virtual contribution. The (j)^ integration is not straightfor- 
ward because the constraint E^r^^^^ < pta does not lead to a simple expression of PrsC'/'s)) 
but when Ei^^^^ is close to px'y we get 

^NLO ^ ^Born ^ ( _ 1 ^^2 Pll_^ I . (1^ 



PTl - ^TM^ 



When -Ejfj^-jj > pt'^ we obtain 



^NLO ^ ^Born ^ I i 1^2 _jEl | . (12) 



2"" M'* 



T,min ~ "PT-i , 



The cross section displayed in fig. § is integrated over pj-^ from a lower limit p^j^jn to 
VT^i = \/s- Therefore the soft NLO contribution to the cross section for E^:^^^ < px-y is 
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where we used the approximation 0"^°'^'^ ps a {pr-y) "'"■'' ^ {''^eff ~ 3). When Ej^^:^^ increases 
up to Pxram^ ^^^ negative contribution ([l3|) to the cross section becomes more and more 
dominant and explains the behaviour of the cross section in the range E^r^^:^^ < Pj-min- 
When -E^'^j^in > Prmin-' ^^^ integral over pTy splits into two parts 

dPT-y = / dpT-y + / . , dpT-f , 

pZ, . JtjL J Ei? ■ 

^r,min ^T,min r,min 



the second part being calculated already in ([131) . The first integral is a positive contribution 
to the NLO cross section, given by an expression similar to (|13|): 
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This positive contribution explains the behaviour of the cross section for Ej, ^^^ > pj, ^ 
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